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Abstract 

We study numerically a dynamic scaling of a diffusion process involv- 
ing the Langmuir type adsorption. We find dynamic scaling functions in 
one and two dimensions and compare them with direct numerical simula- 
tions. We study further the dynamic scaling law on fractal surfaces. The 
adsorption-diffusion process obeys the fracton dynamics on the fractal 
surfaces. 



1 Introduction 

Diffusion and adsorption are fundamental processes in various nonequilibrium 
phenomena, especially in the research fields of fractals PP . The diffusion-limited- 
aggregation (DLA) is a typical fractal pattern generated by the diffusion and 
deposition [2 El ■ The fractal growth and anomalous scaling have been studied 
even in atomistic processes and nanoscale depositions of thin films ^ El E]- 
Transport phenomena on various fractals were also intensively studied, and the 
anomalous dynamical behavior is characterized by the fracton dimension |S] . 
Irreversible reaction-diffusion systems such as A+A-^ O or A+B^ O were also 
studied on fractal substrates E| . 

Adsorption phenomena on fractal surfaces are important in chemistry 
The adsorption on porous media such as activated carbons, zeolites and silica- 
gels is used for elimination, separation, and storage of various materials. The 
fractal dimensions of such porous media were measured by the adsorption of 
particles with various sizes, the small-angle x-ray scattering and the scanning 
tunneling microscopy The transport phenomena including the diffu- 

sion process are important for the temporal evolution of the adsorption onto 
the porous materials. The diffusion and adsorption processes on fractals are 
important also in biology "14". The oxygen is transported to the whole body 
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through ramifying structures of lungs and blood vessels and then adsorption and 
chemical reactions occur at the ends. The fractal dimensions were estimated for 
lungs and blood vessels P El ■ The convective transport by fluid flow plays 
an important role on a larger scale, but the diffusion processes are important in 
smaller scales, because the fluid velocity decays by viscosity. 

In this paper, we study a simple adsorption and diffusion process on regular 
lattices and typical fractals. We consider a reversible adsorption process and 
the total mass is locally conserved in the time evolution. Our model is a dif- 
ferent type problem from irrevesible adsorption or the reaction-diffusion models 
studied before El El ■ We use a deterministic model and do not consider 
statistical fluctuations. We obtain dynamic scaling functions with numerical 
simulations. Although the model and the simulations are very simple, such 
scaling functions for the adsorption process were not reported before. 



2 Model equation and dynamic scaling on reg- 
ular lattices 

Firstly, we consider a one-dimensional adsorption-diffusion model. We consider 
the Langmuir type adsorption |17| . The rate of adsorption is proportional to the 
density c of the diffusion particle and the number of vacant sites for adsorption. 
The occupation fraction 6 {0 < 6 < 1) oi the surface site obeys 

^ = kacil ^6)- kaO, (1) 

where ka and kd are the rate of adsorption and desorption respectively. At 
equilibrium, the fraction 9 is given by — bc/(bc+ 1), where b — ka/kd- The 
diffusion particle obeys a reaction-diffusion equation 

ft=-kac{l-9) + kde + D^, (2) 

where D is the diffusion coefficient. The conservation law oi c + 9 is locally 
satisfied in Eqs. (1) and (2). We study a simple adsorption-diffusion process 
where the particle density c is fixed to a constant value cq at a; = 0, that is, 
the system is in contact with a particle reservoir. The initial conditions are 
assumed to be c{x, 0) = and 9{x, 0) = 0, that is, there are no adsorbates and 
diffusion particles initially. We have performed a direct numerical simulation 
of eqs. (1) and (2) for cq = 1,^^ = l,kd = 0.5, L» = 1 and L 1000. Figure 
1(a) displays three snapshots of c[x,t) at t = 50,200 and 800. The particles 
penetrate into the adsorbent as time goes. Figure (lb) displays five snapshot 
profiles of c{x,t) at t = 200,400,600,800 and 1000 by rescaling the horizontal 
axis by x' = xji^l"^ . The perfect overlapping implies that the time evolution 
of c{x,t) obeys a scaling law c{x,t) = c{x/t^^'^) and 9{x,t) — 9{x/t^/'^). If the 
scaling law is assumed, eqs. (1) and (2) are rewritten as 

, dc „ d^c , , d9 
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where x' = xji^l"^. Substitution of the equilibrium condition Q = 6c/(l + be) 
into eq. (3) yields 

(Pc b \ dc 

- -V2a^ + (1 + 6c)2 ) d^- 

A solution to the nonlinear differential equation (4) approaches for x' — > oo, 



D 



(4) 




Figure 1: (a) Three snapshot profile of c{x,t) a,t t = 50,200 and 800. (b) Five 
profiles of c{x/t^/'^) at t = 200n (n = 1, • • • ,5) (SoUd curves), and a solution 
of eq. (2) (dashed line). They are completely overlapped, (c) Three dynamic 
scaling functions for b = 0.01, 0.5 and 50. 




Figure 2: Seven profiles of c(x/i^/^) at i = 10+50n (n = 1, • • • , 7) (Solid curves), 
and a solution of eq. (3) (dashed line). They are well overlapped. 

only if the derivative dc/dx' takes a special value at x' = 0. The solution can 
be numerically obtained by searching the special value. The dashed curve in 
Fig. (lb) is the special solution, and the solution represents the scaling function 
c(.T/t^/^). If ka = kd = 0, the adsorption and desorption do not occur, then the 
scaling function is exactly solved as 



c{x') 



Co 



f-OO 

/ exp{-x'^/{4D))dx. 

Jx' 



(5) 



The numerically obtained scaling function for ka = I and kd = 0.5 is definitely 
smaller than the scaling function (5), because the adsorption process decreases 
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the density of the diffusion particle. Figure 1(c) displays the dynamic scaling 
functions at 6 = 0.01, 0.5 and 50. As b is increased, the dynamic scaling function 
has a longer tail. The dynamic scaling function at 6 = 50 is very close to the 
function of eq. (5) . 

We have found a similar dynamic scaling law in the two-dimensional system. 
The particle density is fixed to be constant cq at the central point {x,y) = 
{L/2, L/2). We have performed numerical simulations on a square lattice for 
ka = l,kd — 0.5, D = 1 and cq = 1, and the system size is L x L = 600 x 600. 
In the two dimensional system, the last diffusion term in eq. (2) is replaced by 
c{i + 1, j) + c{i ~ 1, j) + c(i, j + 1) + c(i, j — 1) — 4c(j,j), where (i, j) denotes 
a lattice site in the square lattice. Figure 2 displays seven snapshot profiles of 
c{x') at i 10 + 50?! (n = 1, • • • , 7), where x' = [i - L/2)/t^/'^ and j is fixed 
to be L/2. The dynamic scaling law c{i,j,t) — c{r/t^^^) is satisfied for large 
t, where r = — L/2Y + (j — L/2Y is the distance from the center. The 
two-dimensional scaling function is expected to obey 

where x' — r/t^^'^, if the two-dimensional Laplacian is replaced by /dr^ + 
{l/r)d/dr. A numerically obtained solution which approaches to for r — > oo 
is displayed by a dashed line in Fig. 2. The derivative dc/dx' diverges as 1/x' 
at a;' = 0. The theoretically estimated dynamic scaling function is a good 
approximation to the scaling function obtained by direct numerical simulations. 

3 Dynamic scaling of a diffusion-adsorption pro- 
cess on fractal substrates 

We investigate further dynamic scaling laws in the adsorption-diffusion process 
on fractal substrates. The first example is a deterministic fractal set of fractal 
dimension Df = log 5/ log 3 ^ 1.46 as shown in Fig. 3(a), which is called Vicsek's 
fractals JHl- The size of the fractal set we used is L = 3^ = 729, that is, the 
contraction map is repeated 6 times. The diffusion process occurs only on the 
fractal set. Namely, the diffusion term X^i' j'{c(*'i j') ~ c(z,j)} is calculated 
only for the nearest-neighbor sites on the fractal set. Figure 3(b) displays the 
time evolution of the total quantity Q of adsorbed particles The 
quantity grows as where a is close to log 5/ (log 5 + log 3). The exponent is 
evaluated from the relation Q - R^f - (ti/rf™)^/ ^ t^t/d^^ where R is the 
typical radius of the extent of the diffusion particle and d^ is the anomalous- 
diffusion exponent, and the relation R^ ^ f^/du, jg used. The exponent dyj is 2 in 
one and two dimensions. The exponent a is equal to Ds/2, where Dg = 2Df /dn, 
is the fracton dimension. For the loopless fractal, d^ is given by -I- 1 [7lll9|. 
therefore a is evaluated as a = Df/[Df + 1) = log5/(log5-|-log3). That is, the 
dynamic exponent a is determined by the anomalous diffusion on the fractal 
set. This is due to the conservation law for the diffusion particles and the 
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Figure 3: (a) Vicsek's fractals with fractal dimension log 5/ log 3. (b) Time 
evolution of the total adsorbent quantity Q{t). The dashed line denotes the 
curve with exponent a = log5/(log3 + log5) (c) Dynamic scaling of c{x/t^) 
(solid lines) where [3 = l/d^ = log 3/ (log 3 + log 5). The dashed line is an 
approximate function using eq. (7). 

adsorbed particles. However, the form of the dynamic scaling function depends 
on the detail of the adsorption process. Figure 3(c) displays four snapshots of 
cixjt^l'^^") at t = 20000 + 500n (n = 1, • • • , 4), which shows the dynamic scaling 
with exponent (i = 1/dw = log 3/ (log 3 + log 5). If the scahng function of the 
form c{x/t^/'^'^) is assumed, and the diffusion constant D is assumed to behave 
effectively as D(t) ~ D^t^f^-^ IZHI], eq. (4) is rewritten as 

The dashed hne is a solution of eq. (7) for l/d^ — log3/(log3 + log 5) and 
— 1. It is close to the dynamic scaling function obtained by direct numerical 
simulations, but small deviation is observed, probably because the diffusion 
term on the fractal set is approximated by the second derivative D{t)d^c/dx'^ 
with the time dependent diffusion constant. 

The second example is another deterministic loopless fractal of dimension 
2 as shown in Fig. 4(a). The dimension is exactly 2, but there are no loops, 
in contrast to the normal square lattice. Therefore, the anomalous diffusion 
exponent is (i„ = 3, which is definitely different from 2 for the normal two- 
dimensional system, and the fracton dimension is 4/3. Figure 4(b) displays a 
time evolution of the total quantity Q{t) of the adsorbed particles, which grows 
approximately Q{t) ~ t^^^. Figure 4(c) displays four snapshots of c[x/t^/'^"') 
at i = 20000 + 500n (n = 1 , • • • , 4) , which shows the dynamic scaling function 
with exponent 1/dw — 1/3. The dashed line is a solution which satisfies eq. (7) 
for 1/dw = 1/3 and Dq = 1. There is some difference between the solution to 
eq. (7) and the numerically obtained scaling function in this case. 

The third example is an adsorption-diffusion process on a critical percolation 
cluster. The critical site percolation cluster is constructed on a square lattice. 
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Figure 4: (a) Loopless fractal with fractal dimension 2. (b) Time evolution of 
the total adsorbent quantity Q{t). The dashed line is a curve with exponent of 
a = 2/3. (c) Dynamic scaling of c{x/t^) (solid hnes) where [3 — 1/dw — 1/3. 
The dashed line is an approximate function using Eq. (7). 




Figure 5: (a) Adsorbed sites satisfying 9 > 0.01 at t = 100000. (b) Time 
evolution of the total quantity Q{t) of the adsorbed particles. The dashed 
line is a curve with exponent of a = 0.66. (c) Dynamic scaling of c{r/t^) at 
t = 12000,40000 and 100000, where (3 = 1/2.87 



The diffusion and the adsorption occur only on the percolation cluster. The 
density of c is fixed to be 1 at a certain site (200,201). Figure 5(a) displays 
the sites where the adsorbate quantity 9 is larger than 0.01 at t ~ 100000. 
Figure 5(b) displays a time evolution of the total quantity Q[t) of the adsorbed 
particles, which is approximated by Q{t) ~ t^-^^. The fractal dimension of the 
two dimensional percolation cluster is Df = 91/48 and the anomalous diffusion 
exponent is known to be ~ 2.87 ± 0.02 [201 ■ The fracton dimension dg = 
2Df/dw ^ 1.31 ± 0.02, which is close to the Alexander and Orbach conjecture 
ds — 4:/i for the percolation clusters in any dimensions |21j . Our adsorption- 
diffusion process obeys the dynamic scaling law with the same exponent as the 
fracton dynamics on the percolation cluster. Figure 5(c) displays three profiles 
of c(r/ti/2.87) ^ ^ 12000, 40000 and 100000, where r is the distance from the 
center (200, 201) and the average of c is taken for all the sites on the percolation 
cluster, which are located between two circles of the radius r and r + 1. The 
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three plots are overlapping fairly well, which implies that the dynamic scaling 
law is satisfied. 




Figure 6: (a) Time evolution of the total quantity Q{t) of the adsorbed particles 
for the looplcss fractal shown in Fig. 4(a). On the lower-left side, the density 
c is fixed to be 1. The dashed line is a curve with exponent of a = 1/3. (c) 
Dynamic scaling of c{x/t''^) at t = 20000, 25000, 30000, 35000 and 40000, where 
(3=1/3 

We have performed several simulations in which the particle density is fixed 
at a point for the sake of simplicity. However, adsorptive particles usually 
penetrate from the surface of the porous media. We have performed a few 
numerical simulations imposing the fixed value cq on a line. As an example, we 
have used a fractal substrate shown in Fig. 4(a), and the particle density is fixed 
to be c{x, y) = cq at the lower-left side satisfying x+y = 257. On the other three 
boundary sides, no-flux boundary conditions are imposed. We have calculated 
the total quantity Q of the adsorbed particles and the time evolution of c(a;, y) 
on the line satisfying y = 256. Figure 6(a) displays a time evolution of Q (solid 
line), which obeys Q{t) ~ t^^^. The exponent a is reduced by half in contrast 
to the value shown in Fig. 4(b), since the adsorption- diffusion process evolves 
in a quasi-one-dimensional manner. Figure 6(b) displays the dynamic scaling of 
c{x/t^) with the exponent (3 = 1/dw = 1/3, which is the same exponent as the 
one in Fig. 4 (c). The exponent (3 does not depend on the boundary condition. 

We have performed a similar simulation on the random lattice contructed 
by the site percolation used in Fig. 5. The particle density is fixed at the left 
side a; = and the no-flux boundary conditions are imposed on the other three 
sides. Figure 7(a) displays the region with 9 > 0.01 a,t t = 100000, which 
shows the adsorptive particles penetrate into the percolation cluster from the 
left side. Figure 7(b) displays the total quantity Q{t) of the adsorbed particles. 
Q{t) roughly obeys a power law Q{t) ~ with a — 0.33. The exponent a is 
reduced by half, compared to the value shown in Fig. 5(b) owing to a quasi-one- 
dimensional growth from the linear surface. Figure 7(c) displays flve proflles of 
c{x/t^^^-^'^) att = 20000-n (n = 1, 2, • • • , 5), where the average of c is taken over 
the y coordinate on the percolation cluster. The five plots are overlapping fairly 



7 



M (b) (c) 




Figure 7: (a) A snapshot pattern of adsorbed sites satisfying > 0.01 a.t t = 
100000. On the left side a; = 0, the density c is fixed to be 1, in contrast 
to the simulation in Fig. 5. (b) Time evolution of the total quantity Q{t) of 
the adsorbed particles on the percolation cluster. The dashed line denotes a 
line with exponent of 0.33. (c) Dynamic scaling of c{x/t^) at t = 20000 • n 
(n = 1, • • • , 5) with 13 = 1/2.87 

well, which implies that the dynamic scaling law with the (3 — 1/2.87 is satisfied. 
The exponent is the same as the one in Fig. 5(c). Thus, we have found that the 
exponent a of the total adsorption Q{t) depends on the boundary conditions, 
but the exponent fi for the dynamic scaling function does not depend on the 
boundary conditions. 

4 Summary 

We have performed direct numerical simulations of a simple and instructive 
model equation for the diffusion and adsorption process on regular lattices and 
regular two-dimensional fractals, and a two-dimensional percolation cluster. We 
have found the dynamic scaling laws in two kinds of simulations with different 
fixed boundary conditions. The total quantity of the adsorbed particles obeys 
a power law Q{t) ^ t", where the exponent is d/2 for one and two dimensional 
systems and a ^ Df/d^ on fractal sets, when the particle density is fixed at 
a certain point. The number of surface sites (the surface area) is characterized 
by the fractal dimension Dj. It is favorable for the adsorption that the surface 
area is enlarged owing to the fractal surface. However, the diffusion occurs more 
slowly on fractal surfaces than on the smooth surfaces, which is unfavorable for 
the fast adsorption. It is characterized by the anomalous dynamic exponent 
1 /dw The exponent a, which determines the time evolution of the total quantity 
of the adsorption, is given by the product of Df and l/rf,„, when the particle 
density is fixed at a certain point. The exponent a is halved when the particle 
density is fixed at a certain line. The dynamic scaling functions c{x/t^/'^'" ) were 
found on one and two dimensional lattices and several fractals. We have found 
the form of the dynamic scaling function for one and two dimensional lattices 
by solving numerically Eqs. (4) and (5). 
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We have tried to get an approximate dynamic scaling function for the adsorption- 
diffusion process on fractals, however, it is left to future to obtain a better form 
of the dynamic scaling function for the adsorption-diffusion processes on frac- 
tals. Our model is a very simple deterministic model which corresponds to the 
mean-field approximation, however, the applicability to real systems might be 
limited, since statistical fluctuations, interactions between adsorbed particles 
and convective flows are neglected. It is desirable to generalized the simplest 
model to include such various effects. It is also left to study the diffusion- 
adsorption on realistic three-dimensional random fractals using other type of 
adsorption process such as the BET type adsorption |22| . 
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